Quantum periodic cluster methods for strongly correlated electron systems 
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Quantum periodic cluster methods for strongly correlated electron systems are reformulated and 
developed. The reformulation and development are based on a canonical transformation which 
periodizes the fermions in the cluster space. The dynamical cluster approximation and the cellular 
dynamical mean field theory are related each other through the canonical transformation. A cluster 
perturbation theory with periodic boundary conditions is developed. It is found that the periodic 
cluster perturbation theory converges rapidly with corrections 0{1/L'^), where Lc is the linear size of 
the clusters, whereas the ordinary cluster perturbation theory converges with corrections 0{1/Lc). 
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I. INTRODUCTION 



One of the most active areas in the condensed matter 
physics is the search for new methods which could cap- 
ture essential features of electron correlations and disor- 
der in a controllable manner. Perhaps the most success- 
ful and widely used methods are the dynamical mean 
field theory (DMFT)^'- and the coherent potential ap- 
proximation (CPA)i2i^ Both these methods are exact 
in infinite dimensions. However, at finite dimensions 
they neglect nonlocal correlations. This limitation has 
raised interests in developing of quantum cluster meth- 
ods which could capture nonlocal correlations. The most 
successful of self consistent quantum cluster methods are 
the molecular coherent potential approximation (MCPA) 
for disorder systemsj^ and the cellular dynamical mean 
field theory (CDMFT)J and the dynamical cluster ap- 
proximation (DCA)2iS for correlated electron systems. In 
these quantum cluster theories the system lattice is split 
into a series of clusters. Both local and nonlocal corre- 
lations within the clusters are treated exactly, whereas 
the nonlocal correlations between different clusters are 
treated in a mean-field approximation. The mean field is 
taken into account by mapping the lattice problem onto a 
self-consistent effective cluster problem. The MCPA and 
CDMFT are formulated on the real space, and share a 
common microscopic definition. Within the CDMFT (or 
the MCPA) the cluster Green function is calculated with 
open boundary conditions. The DCA is traditionally for- 
mulated in the reciprocal space. It is based on the idea of 
discretizing irreducible quantities on the reciprocal space. 
Within the DCA the cluster Green function is calculated 
with periodic boundary conditions. Due to the difference 
of boundary conditions it seems that the DCA and the 
CDMFT (or the MCPA) have different microscopic def- 
initions. For instance, the CDMFT can be formulated 
within the self-energy functional approach, whereas it 
seems that the DCA cannot beii2*ii There is also a view 
of the DCA in the real space representation which shows 
a relation between the DCA and CDMFT The 
view is based on an unitary transformation of the hop- 
ping matrix. However, the interaction part still remains 



unchanged. In this paper we show that the DCA and the 
CDMFT share a common microscopic definition through 
a canonical transformation. The canonical transforma- 
tion is similar to the unitary transformation which was 
used to view the DCA in the real space. It periodizes 
the fermions in the cluster space, that the cluster Green 
function has periodic boundary conditions. We apply the 
canonical transformation to the whole Hamiltonian, and 
in the presence of a constrain which prevents the umklapp 
momentum transfer from the superlattice to the cluster 
space we obtain a periodized Hamiltonian. The umk- 
lapp momentum transfer appears when the sum of two 
momenta of the superlattice is beyond the first Brillouin 
zone, and is solely due to the imaginary superlattice con- 
struction of clusters with periodic boundary conditions. 
Within the periodized Hamiltonian the DCA can be de- 
rived from the CDMFT. In the such way, the DCA and 
the CDMFT can be unified into a common microscopic 
background. The derivation also gives the microscopic 
background of the discretization of irreducible quantities 
on the reciprocal space, and clarifies the approximation 
nature of the DCA. 

The present paper consists of two parts. In the first 
part a derivation of the DCA from the CDMFT through 
the canonical transformation for fermions is presented. 
The second part is concerned with the cluster pertur- 
bation theory (CPT),iSiiSJi The CPT is also a quan- 
tum cluster approach. However, in difference with the 
CDMFT or the DCA, the CPT does not have the self 
consistency. The CPT can be viewed as the first term 
of a systematic expansion around strong coupling. How- 
ever, within the CPT the cluster Green function is cal- 
culated with open boundary conditions, that the wave 
vector within a cluster is no longer a conserved quantum 
number. There are also several approaches which gener- 
ate periodic boundary conditions for the CPT. One is to 
add appropriate hopping terms and then subtract them 
within strong coupling perturbation theory.^* However, 
this approach give results which are less accurate than 
the ones from open boundary conditions.^'' The other 
approach adapts the periodicity of the clusters from the 
DCA. It was used to study the Ising modeliS In this pa- 
per we present the periodic cluster perturbation theory 
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for fermion systems. The periodic cluster perturbation 
theory (PCPT) is the CPT formulated for the periodized 
fermions which are obtained after making the canonical 
transformation. The PCPT has truly the small parame- 
ter 1/Lc, where Lc is the cluster linear size. It turns out 
that the PCPT converges quadratically in 1/ Lc, whereas 
the CPT converges linearly in l/Lc- This is also tested 
on the study of an exact solvable model. 

The plan of the present paper is as follows. In Sec. II 
we present the canonical transformation which trans- 
forms the original fermions onto periodized ones. In 
Sec. Ill we derive the DCA from the CDMFT after mak- 
ing the canonical transformation. The PCPT is pre- 
sented in Sec. IV. In this section we also study an ex- 
act solvable model in both the direct and periodized 
forms. Finally, the conclusion and remarks are presented 
in Sec. V. 



II. PERIODIZED FERMIONS 

We present a periodization for fermions through a 
canonical transformation of fermionic operators. For an 
illustration a fairly general model of correlated electron 
systems is considered. The Hamiltonian of the model 
reads 

where {a} = {cti, (72, cs, 0-4}. c^^ (Cj^.) is the creation 
(annihilation) operator for fermion with spin a at lattice 
site i. t{i,j) is the hopping integral, and U^'^^{i,j) is the 
interaction strength. We will consider a hypercubic lat- 
tice of linear size L with the lattice constant a = 1 on d 
dimensions. The number of lattice sites thus is N — L"^. 
The lattice is divided into a set of identical clusters of lin- 
ear size Lc- The number of lattice sites inside a cluster 
thus is Nc = L'^. The set of the clusters form a super lat- 
tice. We use letters i, j,--- to label the lattice site of the 
original lattice, and use to denote the position of site 
i. Letters k, p, ... are used to denote the wave vectors of 
the reciprocal space of the original lattice. Capital letters 
/, J,... are used to label the origin of the clusters. The 
origin coordinate of clusters are denoted by R/, and the 
wave vectors of the reciprocal space of the superlattice are 
denoted by K, P, ... The lattice sites inside a cluster are 
labeled by a, 6,... Their position is denoted by Fa, Ff,,... 
and the wave vectors of the reciprocal space of the cluster 
are denoted by k, p,... Thus = R/ -I- fa, k = K -I- k,... 
Note that exp(ik • R) ~ 1, since ka — StTOq/Lc, and 
Ra = niaLc, where ria, rua are integer. We use the 
Greek letter a to denote the spin (or other quantum) 
variables. With these notations we rewrite the operators 
and model parameters as Cia = Caa{I), t{i,j) = tab{I, </), 
C/{'^>(i,j) = C/];^(/, J),... Denoting 4^, c^,, the Fourier 



transform of cj^ , c^^ , the Hamiltonian |^ can be rewrit- 
ten in the reciprocal space representation 

i/ = ^t(k)cLck. 

k,(T 

+^ E Ef^^'^^(p)4+p,..Ck,.,4_p,.3Ck,.„(2) 

k,k',p {a} 

where i(k) and f7^°'^(p) are the Fourier transformation 
of t(i,j) and {/^"^{ijj), respectively, i.e., 

^(^'^■) = ^E^(k)e^''-('--r.), (3) 

k 

C/M(z,j) = l^[/W(k)e'''-('-.-r.). (4) 

k 

We introduce operators 

C-a.(K) = yK^,^^(J)e-K.(R.+.J^ 

c-L(K) = y^^cL(/)e^K,R.+..), (g) 
/ 

One can check that these operators satisfy the anticom- 
mutation relations 

{Ca.(K),4,(P)} = SabS„,,d{K-P), 
{Ca<.(K),Cfc,,(P)} = 0, 

{c-L(K),4,(P)} = 0, 

where {A, B} — AB + BA. Since the anticommuta- 
tion relations are preserved, the transformation ®-® 
is canonical. The inverse transformation reads 

ca.{l) = y^5:c-a.(K)e^K-(^^+^»), (7) 

K 

cUl) - Ec"L(K)e-^^-(^^+^"^ (8) 

K 

We also introduce other operators in the real superlattice 
space 

-CaM) = \/¥E^-(K)e^^•^^ (9) 

K 

cUl) - E^1^(K)e-'^-^^ (10) 

K 

One can check that these operators (P)l- l(Tn|) also satisfy 
the fermionic anticommutation relations. They annihi- 
late or create a fermion at site (/, a) . The transformation 
(I^- (I1U|I can be viewed as the Fourier transformation for 
the fermion operators Caa (K) , (K) in the superlattice 



3 



space. From Eqs. (|^- (|10() and jSJ-lEll, one can show that 

-CaM) - /^E^K+k,.e^''-''^+^^--, (11) 

--UD = /^E4.k.^-^"-"^-"-"- (12) 

K,k 

In deriving Eqs. pi|) - (|12|l we have used exp(ik-R/) = 1. 
Since exp(ik • fa+L^) = exp(ik • fa), one can see that 
the fermionic operators Caa {!) , (/) are periodic in the 
indexes a, whereas the original ones c^^{I), c\^{I) are 
not. The fermion created by c\^{I) is periodic in both 
the superlattice and cluster spaces. Due to the periodic 
property we call Caail), the periodized fermion 

operators. The periodicity allows to perform the Fourier 
transformation of Caa{I), in the superlattice and 

cluster spaces 

C£,(K) = yiEc..(/)e-'^K-^-'^-^«, (13) 

I ,a 

4(K) = J^Y.^Ul)e''^-'''+'^-''- (14) 

I,a 

We want to emphasise here that both K and k should 
be restricted to the first Brillouin zones of the superlat- 
tice and cluster spaces. From Eqs. (Ill|) - (|12|l we obtain 
immediately c^JK) = c^^^^^^, c^JK) = ct^j^_^. How- 
ever, these relations are valid only for momenta in the 
first Brillouin zones. Thus, from the Hamiltonian ||2Jl we 
obtain 

H^Y. i(K + k)4(K)c-,(K) 

K,k,(T 

+^ E E Ef^^^^(p+p)4+p,..(K+p) 

K,K',Pk,k'^p {a} 

cr (K)ct _ (K'-P)cj-, (K'). (15) 

Note that in Hamiltonian H15() all sums over momenta 
are restricted to the first Brillouin zones. Moreover, 
when K + P and K' — P are beyond the first Bril- 
louin zone, one has to translate them back to the first 
Brillouin zone, thus the momentum sums still are over 
the full first Brillouin zone. Equivalently, the momenta 
of the periodized fermion operators take values Ka = 
27rmod(nQ,, L/Lc)/L — tt/Lc. In deriving the Hamilto- 
nian (|15|) we have used the following restriction 

(5(p + k-q) = 5(P + K-Q + G)(^(p + k-q + G') 
= ,5(P + K-Q + G),5(p + k-q), (16) 

where G, G' are arbitrary momenta which are either 
zero or a vector of the reciprocal superlattice space. G, 
G' are also the momenta of the reciprocal cluster space. 



G 7^ just means an umklapp process on the superlat- 
tice, and G' can be interpreted as a momentum trans- 
fered by the umklapp process from the superlattice to 
the cluster space. The first line of Eq. (|16|) is exact, 
while the second line is a particular restriction. This re- 
striction prevents the momentum transfer in the umklapp 
processes of the superlattice. The umklapp processes ap- 
pear when the sum of two momenta of the superlattice is 
beyond the first Brillouin zone, and they transfer a mo- 
mentum modulo a vector of the reciprocal superlattice to 
the reciprocal cluster space that generate a momentum 
nonconservation. However, the such momentum transfer 
is unphysical, it is solely due to the imaginary discretiza- 
tion of the original lattice into superlattice of clusters 
with periodic boundary conditions. In order to maintain 
the correct physics, one has to impose a constrain which 
forbids the umklapp momentum transfer on the super- 
lattice. The restriction (|16|l is the constrain of the super- 
lattice construction. As we will see later on an example 
of exact solvable model, the constrain indeed leads the 
periodized Hamiltonian H15|l to the exact solution. Per- 
forming the Fourier transformation ()13|I - H14|I . from the 
Hamiltonian H15|) we obtain 

E u!^:Hl,J)cLAI)^aaM)cU{m.AJ), (17) 
IJ,a,b,{a} 

where the hopping integral and interaction of the peri- 
odized fermions now are 

K,k 

ulbHl, J) = ^J2 U^'Hi^ + k)e'K-(i^^-i^^)+^k.(r.-f,)_ 

K,k 

(19) 

One can notice that iab{I , J), uj^^\l, J) are cyclic in the 

indexes a, 6, whereas tab{I,J), U^l^{I,J) are not. This 
means that the Fourier transformation of tabiI,J) and 
uj^b^ (I, J) are diagonal in the reciprocal space. Denoting 
f£(K), U^'^\k.) the Fourier transformation of tabil, J), 
uj^b^ {I, J), one can see immediately from Eqs. (|18ll - (|19|) 
that tk(K) = <(K + k), ul'"^K) = [/{^>(K + k). Note 
that the hopping term H18(l was also obtained previously 
by employing an unitary transformation. ^^'^"^ The uni- 
tary transformation is similar to the canonical transfor- 
mation for fermion operators. Beside the periodicity, the 
periodized hopping p8|l and interaction H19|l are quite 
different from the original hopping and interaction, re- 
spectively. They may connect those lattice sites that the 
original ones do not. For instance, if the original hop- 
ping is nearest neighbor then the hopping only couples 
the nearest neighbor sites either within the clusters, or 
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on the cluster boundaries of the nearest neighbor clus- 
ters. However, in the same case the periodized hopping 
can couple a lattice site with any other lattice site. The 
Hamiltonian (|17l) is adequate to quantum cluster meth- 
ods with periodic boundary conditions. 

Next, we consider the one-particle Green function of 
the periodized fermions 



cluster is 



Ga6(K,z) =« Ca.(K)|4(K) >3 . 

Without difficulty one can show that 

Gab{K) = ^ ^ G(K + k)e^i^ (^"-^^), 



(20) 



(21) 



where G(k) =<C Ckcrjcj^.^ 3>2 is the Green function of the 
original fermions. Equation H21|) shows that the Green 
function of the periodized fermions is diagonal in the re- 
ciprocal cluster space. We immediately obtain 



G(K + k,^) = Gk(K,z), 



(22) 



Equation l|22|) expresses the identity of the Green func- 
tions of the original fermions and of the periodized 
fermions. The Dyson equation of the Green function 
Gat(K, z) reads 



[G(K,z) 



zi -t{K) -t{K,z) 



(23) 



where the hat symbol denotes the matrix notation in the 
cluster space. S(K, z) is the self energy. Since the Green 
function G(K, z) and the hopping f(K) are diagonal in 
the reciprocal cluster space, the self energy S(K, z) must 
diagonal in the reciprocal cluster space too. This is con- 
sistent with the Hamiltonian H17() . where both the hop- 
ping and interaction terms are diagonal. 



III. DYNAMICAL CLUSTER 
APPROXIMATION 

With the periodized Hamiltonian (|17|l we can derive 
the DCA from the CDMFT (or from the MCPA). Ap- 
plying the CDMFT to the periodized fermions we obtain 



G(K,z) = [zi ~i{K) ~j:{z)\ 



(24) 



The hopping matrix i(K) is the Fourier transformation 
of the periodized hopping (|18() in the superlattice space, 
i.e., 



U(K) = ±Y.t{K + k)e^- 



N, 



(25) 



The self energy E(z) is determined from an effective sin- 
gle cluster problem. The action of the effective single 



^eff = - r f dTdT'Y,cUT)g-,\T~T')c,AT') 

+ Ull\ci,^-CaaAa,Cb^.){r), (26) 

a, 6, {a-} 

where Q^]^{t — r') plays the role of the effective mean 
field acting on the cluster. The interaction in the effective 
single cluster is obtained from Eq. I|19|) 



K,k 



(fa-rt) 



(27) 



The self consistency requires the identity of the cluster 
Green function obtained from the effective single clus- 
ter problem (|26|l and the cluster Green function of the 
original cluster, i.e.. 



N 



^G(K,.)=[r {z)-t{z)] 



(28) 



K 



Since in the reciprocal cluster space the lattice Green 
function G(K,z) and the self energy S(z) are diagonal, 

the Green function of the effective medium Q {z) must be 
diagonal too. After making the Fourier transformation, 
we obtain from Eqs. (|23I), it^ 



Gk(K,z) 



z-fk(K)~I]k(z)' 



A^4- ' ^ g^\z)-^-^{z) 



(29) 
(30) 



S_ince t^[K) = t_(k+K), t/g = [NJN) C^(K+k), and 
Gj^(K, z) — G(k -I- K, z) we recover the DCA equations. 
Equations H24() - (|28|l can be viewed as the DCA equa- 
tions formulated in the real cluster space. Originally, the 
DCA is formulated in the reciprocal spacei^i It is based 
on the idea of discretizing irreducible quantities on the 
reciprocal space. Later, the DCA was also viewed in 
the real space by introducing an unitary transformation 
for the hopping matrixii2iiii4 However, in this view only 
the hopping term is transformed, whereas the interaction 
term remains unchanged. As a consequence, the hopping 
and interaction terms are not treated on an equal footing. 
It correctly recovers the DCA only for a local single-site 
interaction. In the presented approach both the hopping 
and interaction terms are treated on an equal footing. 
It also clarifies why within the DCA the lattice quanti- 
ties, for instance the hopping term in the lattice Green 
function (Eq. H25I) '). are unchanged, whereas the cluster 
quantities, for instance the cluster interactions (Eq. l|77|) '). 
are discretized or coarse grained. Thus, the DCA can be 
used in both the reciprocal and real spaces. In the real 
space the DCA is just the CDMFT applied to the pe- 
riodized Hamiltonian (|17|) . Due to this formulation the 
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DCA has the properties of the CDMFT. For instance, 
the DCA can also be formulated within the self-energy 
functional approach in the real space by the same way as 
of the CDMFT .iaiii So far we have shown that both the 
DCA and CDMFT share a common microscopic defini- 
tion. They are related each other through the canonical 
periodization transformation. An extensive discussion 
about the two methods is already given in the recent 
review^ 



IV. PERIODIC CLUSTER PERTURBATION 
THEORY 

The CPT is also a quantum cluster approach, although 
it does not have the self consistencyii^*iSiii In this sec- 
tion we present a CPT with periodic boundary condi- 
tions. We refer it as PCPT. Basically, the PCPT is the 
CPT applied to the periodized fermions. A similar PCPT 
which employs the formulation of the DCA for spin sys- 
tems was used previously to study the Ising modeli^ We 
will consider the system with the nearest neighbor hop- 
ping. In the periodization version the hopping H18|) is 
no longer only nearest neighbor. It may couple a lattice 
site with any other lattice site. As the CPT, the PCPT 
splits the intra- and inter-clustcr parts of the hopping 
term. The intra-clustcr part of the Hamiltonian is 

H = igbC^^Cba + ^ifc ^jj^.^ Cao-2 ^^0-3 ^6o"4 i (31) 

a.h.a a,6,{cr} 

where U = iab{IJ), and f/^^"^ = ul^i\lJ)- The inter- 
cluster part is 

H'^= J2 Stab{I,J)ciM)cta{J), (32) 
I ,J,a,b,(7 

where Siab{I, J) = tab{I, J)— tab- Here we have neglected 
the inter-cluster interactions. Originally, the CPT is for- 
mulated for local single-site interactions, thus the inter- 
cluster interactions are absent. Often quantum cluster 
approaches either neglect the inter-cluster interactions or 
treat them in a static mean-field approximation. How- 
ever the inter-cluster interactions can be incorporated 
beyond the static mean-field approximation by employ- 
ing a self consistent approach>i2i2S Note that the cluster 
hopping tab may couple any sites within the cluster. At 
a finite cluster it is reduced by a factor of 0{1/L'l)& 
The inter-cluster hopping St ^ ^/Lc for large linear clus- 
ter size Lc compared to its counterpart t of the original 
CPT.^-^ The nature of the approximation of the CPT is 
a strong coupling expansion in which the small param- 
eter is the inter-cluster hoppingji^^iSiii Therefore in the 
PCPT the small parameter is truly l/Lc- 

Within the PCPT the cluster Hamiltonian (jHU is 
solved exactly. We obtain the Dyson equation for the 
cluster Green function 



GUz) 



1 



where = (Nc/N) J^k ^(K-l-k), and S^(z) is the cluster 
self energy. Here the Dyson equation is written in the 
reciprocal space representation. We can make the Fourier 
transformation in the cluster space because the fermions 
are periodized in the cluster space, thus the cluster Green 
function and self energy are diagonal in the reciprocal 
space. As the CPT, within the PCPT the self energy of 
the lattice Green function is approximated by the cluster 
self energy. Thus we obtain 



Gk(K,z) 



1 



tk(K) - S]£(z) 



(34) 



Since Gk(K, z) = G(K + k, z) and t^,{K) = i(K + k) we 
finally obtain the Green function of the original fermions 



G(K-f k,z) 



1 



z-t{K + k) 



(35) 



Equations are the principal equations of the 

PCPT. As the CPT, the PCPT is exact in the hmits 
Lc — > 00, U/t = 0, and t/U — for local interac- 
tions. The PCPT is formulated in the reciprocal space, 
hence it avoids the diagonalizations of the Green func- 
tion in both real and reciprocal spaces, that significantly 
reduces the computational time compared to the CPT. 
Note that the present PCPT is different to the ordinary 
periodic CPT,i2ii§. where the periodicity is taken into 
account by including additional hoppings. In the ordi- 
nary periodic CPT the hopping is still nearest neighbor, 
while in the present PCPT it is no longer only nearest 
neighbor. Moreover, in the present PCPT all hoppings 
within a cluster are treated exactly, whereas in the ordi- 
nary periodic CPT the adding and subtracting hoppings 
are treated on different footings. 

For an illustration we study the one dimensional large- 
Af model, originally introduced by Afffeck and Marston 
in two dimensions. The model was also used for com- 
parison between the DCA and CDMFT schemesiiSiH We 
will compare the PCPT and CPT schemes to its exact 
solution. The Hamiltonian of the model reads 

<i.j>.cr i(T 

+ ^ 4c..'SV'S-' (36) 

<i,j> .a.{T' 

where < i,j > denotes the nearest neighbor sites, and 
(T = 1, ...,J\f. As usually, t is the hopping integral, J is 
the exchange strength, and /i is the chemical potential. 
We will consider the large TV limit. In this limit the 
model can be solved exactly. Indeed, the quantity x — 

(1/7V) X)o- 4o-'^i+i o- does not fluctuate, that the static 
mean field theory is exact. We obtain the exact solution 



(33) 



^ = — 2^/(/3i?(k))7(k), (37) 

k 

E{k) = -(t + Jx'=''-')7(k)+M, (38) 




FIG. 1: The dispersion Ep(P) of the periodized fermions (the 
open symbols) for the cluster size Lc = 4 in comparison with 
the dispersion E(P + p) of the original fermions (the line). 
The filled symbols are i5p(P) of the umklapp case for Lc = 4. 
{t = J = 1, fi = 1, /J ^ 10). 



FIG. 2: The dispersion Ep(P = 0) of the periodized fermions 
for various cluster sizes Lc (the open symbols) in comparison 
with E(P + p) of the original fermions (the line) at P = 0. 
The filled symbols are -^^^(P = 0) of the umklapp case. The 
inset plots x" of the umklapp case in comparison with the 
exact value, (t = J = 1, ^ = 1, /3 = 10). 



where /3 is the inverse temperature, and f{x) = 
l/(exp(a;) + l) is the Fermi function, and 7(k) = 2cos(fc). 

Within the CPT, the cluster of size has the following 
Hamiltonian 



a—l cT 



a—l CT 



2J^ ^ ^ '^l,crCa,cr'cl+i,CT''^a+l,cr + (39) 



a—l a, a' 



The cluster Hamiltonian H39() can be solved exactly, for 
instance by the exact diagonalization. However, at the 
large-A/" limit it can simply be solved exactly. Indeed, 
the cluster quantity xl = (V-^) 4,<TCa+i,<T does 
not fluctuate, and the cluster Hamiltonian H39I) becomes 
quadratic and can simply be diagonalized. However, in 
this case Xa depends on the cluster index. The clus- 
ter Green function is 



G%z) = [{z-f,)i-i-i:T\ 



(40) 



where t is the cluster hopping matrix, and S*^ is the clus- 
ter self energy 



^ab — -Jixl^a,b-1 + Xb^a,b+l)- 



(41) 



The cluster quantity Xa be calculated from the clus- 
ter Green function 

= - / ^/(/3^)ImG^+i,a('^ + *0+). (42) 

Equations (|40|) - (|42|l give the exact solution of the cluster 
Hamiltonian H39|) at the large-A/" limit. Once the exact 



cluster solution is obtained, the lattice Green function 
within the CPT is computed bjiiSJI 



G' 



CPT 



(k,z) - ^^Gab(K,2)e"^''-(^ 



(43) 



ab 



where G~l{K, z) ^ {z - ^i)5ab - U(K) - Sjjj,. The lat- 
tice quantity x = (l/-^^) X^ct ^\cr^i+i,cr within the CPT is 
calculated from this lattice Green function 



CPT 



duj 



-/(/3c.)ImG,%(c. + zO+) 



= -zE/ ^/(MlniG„+i.,(K,c. + zO+).(44) 

K,a 

Next, we consider the one dimensional large-A/" model 
(|36|l in the periodization form. The periodized Hamilto- 
nian reads 

H= tab{I,J)cUl)CbAJ)+t^T.^aAl)CaAl) 

I .J,a,b.{T I .a.(T 

+^ E Jab{i,j)ciMYc,..{i)iAJ)^ba{j)m 

I ,J.a.b,a.{T' 

where tab{lTJ) and Jab{I,J) are determined by the 
periodization transformation (|18|l and H19|l . respec- 
tively. The Hamiltonian (|45l) can be solved exactly 
in the large-TV" limit. Introducing quantity Xp(P) = 
(l/A^)E.4<x(P)Cp.(P) we obtain 

Xp(P) = /(/3^p(P)), (46) 
^p(P) = -t7(P + p)-f/i 

-ji^7(P + p-K-k)xk(K). (47) 



K,k 



When Eqs. lj^ - l|T7|) are solved, the quantity x = 
(1/7V) Ylicr '-lo-'^j+i a ^^"^ original fermions can be cal- 
culated too 



(48) 



P.p 



Numerically, we alway obtain — 1^'='^^'=* independently 
on the cluster size Lc- Moreover, since the Green function 
<Ck(K)|ct(K) >=<C c(K + k)|c'l'(K + k) >, one must 
have the identity of the dispersions of periodized and 
original fermions £'p(P) = E{V + p). In Fig. ^ we plot 
the dispersions for the cluster size = 4. It confirms 
the identity of the dispersions of periodized and original 
fermions. In Fig. [5] we plot the dispersions at fixed P = 
for various cluster sizes Lc- It shows the identity of the 
dispersions of periodized and original fermions is inde- 
pendent on the cluster size. These results confirm the 
equivalent of the periodized and original Hamiltonians. 
Next, we consider the umklapp case where the constrain 
is relaxed by its exact relation. Then we obtain 



Xp(P) 



/(/3^p(P)), 
-t7(P + p) 



-^iE7(p 

K,k 



K 



(49) 



k + G)x^(K), (50) 



where G is the umklapp transfer momentum. Since 
G| ^ 1/ic, in the large Lc limit this umklapp case 
approaches to the normal one. The lattice quantity 
X = (l/A/") X^cr o-^i+i a of the original fermions can be 
calculated by 



(51) 



In numerical calculations we take G — 211 /Lc- It turns 
out numerically that x" is never equal to ^ although 

it approaches to the exact value as l/L^, as shown in 
the inset of Fig. |21 In Fig. and [21 we plot also the 
dispersion i?p(P). These figures show that the dispersion 
of the periodized fermions in the presence of the umklapp 
momentum transfer is not exact. These results confirm 
that the constrain is important in order to maintain the 
correct physics. 

Now we apply the PCPT to the one dimensional large- 
TV model (ll-tfi|) . The cluster periodized Hamiltonian is 
obtained from (|45|l 



+ 



a,b,(7 

J 



2M 



E TabC^C^.x'cL'Ch.T- (52) 



At the large-A/" limit the cluster Hamiltonian (|52|) can 
simply be solved exactly. Introducing the cluster quan- 
tity x| — {^/■^)J2a^t.cr^kcr obtaiu the cluster self 
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FIG. 3: The lattice quantity x as a function of inverse tem- 
perature P for the CPT and PCPT with different cluster sizes 
(t ^ J ^ 1, fi = 1). 



energy 



P 



(53) 



where 7^ = (Lc/L) J2k.^(^ + "^^^ cluster quantity 
X^ can be calculated from the cluster Green function. 
We obtain 



41: 



7k-pXf 



(54) 
(55) 



Finally, once the cluster equations (|54|I - H55I) are solved, 
we obtain the lattice Green function of the original 
fermions 



GPCPT(k,z) = G£(K,z) 



1 



z + t7(k)-M-S5 



(56) 



The lattice quantity x^^^'^ = (1/7V) J2a 4,a(^i+i,a is cal- 
culated from this lattice Green function. We obtain 



X 

pPCPT 



PCPT ^ J_^^(^^PCPT(k))^(k), (57) 

(58) 



k 



So far we have obtained the quantity x in both the CPT 
and PCPT, as well as its exact value. We use it to com- 
pare the error and the convergence of the CPT and PCPT 
schemes. 

In Fig. O we plot the lattice quantity x as a function 
of inverse temperature calculated within the CPT and 
the PCPT. It is clear that both the CPT and PCPT 
results converge well to the exact solution. The CPT 
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results converge consistently, whereas for small cluster 
sizes the PCPT results fluctuate around the exact solu- 
tion. However, the PCPT results converge faster than 
the CPT ones. Even at small cluster sizes the PCPT re- 
sults are closer to the exact solution than the CPT ones. 
In Fig.Qlwe plot x as a function oil/ Lc at fixed temper- 
ature. It shows that the CPT results converge linearly 
in \/Lc^ whereas the PCPT results converge quadrati- 
cally in 1/ Lc- The convergence features of the CPT and 
PCPT are quite similar to the ones of the CDMFT and 
DCAi2ii2ii^ The convergence can be understood in the 
term of the hybridization in the locator expansiounSi In- 
deed, the convergence of x depends mostly on the conver- 
gence of the coarse grained Green function (see Eq. (02)) 



canonical transformation to the Hamiltonian, the peri- 



(59) 



K 



where the Green function G'(K, z) can be rewritten as 
G(K,z)=[[G^(z)]-i-<5£(K)]-\ (60) 

This is also valid for the PCPT just by replacing the 
coarse grained Green function of the original fermions by 
its counterpart of the pcriodizcd fermions. The coarse 
grained Green function can be rewritten in the locator 
expansion as^ 



G(z)^[[G^(z)]-i-f(z)]-\ 
where T[z) is the hybridization function 

L,. n-l 

K 

^^f(K)G(K,z)^f(K). 



(61) 



f(z) = [/+^5]<5i(K)G(K,z) 



L 



(62) 



K 



The hybridization is an amplitude for fermion hopping 
from a cluster into the surrounding and back again into 
the cluster. It acts as the small parameter in the expan- 
sion of the coarse grained Green function (|61|l . As in the 
DCA, the inter-cluster hopping 5t in the PCPT scales 
like 1/Lc for large cluster sizes, whereas in the CPT and 
the CDMFT it is of order 1. The average hybridization 
per cluster site scales like l/i^ the PCPT, and like 
\/Lc in the CPTi^ Therefore the PCPT converges like 
\/Ll, whereas the CPT converges like 1/Lc- Note that 
these convergences are valid for the coarse grained Green 
function in any dimensions. The convergence feature is 
an advantage of the PCPT in comparison with the CPT. 



V. CONCLUSION 

We have introduced a canonical transformation which 
periodizes fermions in the cluster space. By applying the 
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FIG. 4: The lattice quantity x as a function of 1/Lc for the 
CPT and PCPT at fixed temperature /3 = 10 (t = J = 1, 



from the superlattice to the cluster space. Within the pc- 
riodizcd Hamiltonian, the DCA can be derived from the 
CDMFT. In the such way, the DCA and the CDMFT 
can be unified into a single microscopic definition. It 
also gives an alternative microscopic background of the 
discretization of irreducible quantities on the reciprocal 
space, and clarifies the approximation nature of the DCA. 
We also develop the PCPT in which the wave vector 
within a cluster is a conserved quantum number. This al- 
lows to avoid the matrix inversion of the Green function 
that reduces significantly the computation time when the 
cluster size is large. The PCPT can work on both the di- 
rect and reciprocal spaces. As the CPT, the PCPT is 
exact in the limits Lc oo, U/t — 0, and t/U = 0. It 
is also clarified that the small parameter of the PCPT 
is 1/Lc- As a benchmark the exact one dimension 1/JV 
model is studied. It turns out that the PCPT converges 
rapidly with corrections 0{1/L'c), whereas the standard 
CPT converges with corrections 0(1/Lc). 
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